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1. Introduction 

A wide class of problems concerned with the quantum information processing requires an 
efficient and robust methods of the creation of highly-entangled states from individual 
qubit states [1], [2]. Therefore, the search and the development of reliable schemes 
for entanglement production appears to be a rather important quantum computational 
issue. Recently, several proposals for the entanglement generation in the semiconductor 
nanostructures based on the quantum dots (QDs) (see, e. g., Refs. [3] - [6]) and its 
application (see, e. g., Refs. [7j - [13]) have been made. In those works, the fault-tolerant 
quantum algorithms [7j - [9], the mechanisms of the quantum teleportation - |12j . 
and the cluster state preparation in the measurement-based quantum computations [T3] 
exploit specific entangled charge or spin states of electrons bound in the QDs. 

In particular, the logical states of an individual (physical) qubit can be presented by 
two single-electron orbital states localized in the spatially separated potential minima 
of a double quantum dot (DQD) p3] - p2]. The quantum operations on that so- 
called charge qubit can be accomplished via adiabatic variation of the DQD confinement 
potential [H] - [17] or by the electromagnetically induced resonant transitions between 
the DQD states [H] - [22]. In order to perform some non-trivial two-qubit operation 
on an arbitrary pair of qubits in a quantum register one should organize the interaction 
between those qubits for a finite period of time. In general, after interaction is off, 
the qubits become entangled with each other, i. e., their total wave function cannot 
be presented as the product of an individual qubit wave functions. The charge qubit 
entanglement via electrostatic control over the tunnel coupling between two neighboring 
qubits has been studied in Refs. [i]. [6]. pi]. [11] - [T7]. 

In this paper, we suggest another technique of entanglement production in the two- 
dimensional array of the charge qubits. The coupling between the qubits is now indirect 
and mediated by an auxiliary structure (AS) fabricated from the quasilinear chains of 
the QDs and containing a single (probe) electron in the quantized part of its conduction 
band. The single-electron energy spectrum of that structure is affected by qubit states 
through the electrostatic interaction between the probe and qubit electrons. In order 
to demonstrate how to entangle the qubits we analyze in detail the probe electron 
dynamics in a separate quasilinear QD AS capacitively coupled to the charge qubit. As 
we shall see, appropriate variation of the confinement potential of such an AS through 
the application of compensating voltages aligns the energy levels of individual AS QDs 
for one of two logical qubit states. It amounts to sharp increase of the tunnel coupling 
between the AS QDs that, in its turn, allows for efficient probe electron transfer (PET) 
along the AS through one of the hybridized AS states (the transport state). The transfer 
protocol is accomplished via the resonant transitions connecting two states each localized 
in corresponding edge AS QD and the transport state delocahzed over the AS. The 
numerical simulation of electron dynamics under tight-binding approximation confirms 
the possibility of successful PET implementation. 

Two points are important here. Firstly, PET from one edge AS QD to another and 
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back results in the implementation of the phase operation on attached qubit. Secondly, 
the PET between the edge AS QDs takes place only if the qubit is in the predetermined 
logical state. Therefore, the following dynamics of probe electron appears, in its turn, to 
be predetermined (conditioned) by the qubit state. It enables us to organize conditional 
quantum operations (in particular, controlled-phase operation) that is the key part 
of entangled state preparation. Using these results and assuming that PET in a 
complex planar AS can be divided into the sequence of elementary PETs along its 
quasilinear components, it becomes possible to construct any desired entangled state of 
the charge qubits attached to that AS. Particularly, the set of described manipulations 
underlying controlled-phase operations, combined with single qubit Hadamard rotations, 
is sufficient for the nine-qubit Shor encoding procedure [23]. Moreover, the conditional 
PET may be exploited for an error syndrome measurements and for usual qubit state 
measurement as well. 

The paper is organized as follows. In Section 2 we examine in detail the energy 
spectrum of the auxiliary QD structure attached to single qubit as well as the 
organization of conditional probe electron dynamics in that structure. This specific 
operation amounts to a single-qubit phase shift and can be used for implementation 
of controlled-phase shift of target qubit in two-qubit circuits and, therefore, for 
entanglement generation. Section 3 describes the algorithm of nine-qubit Shor encoding 
which exploits two-qubit gates based upon conditional probe electron evolution to 
entangle the qubits. Auxiliary structure used in the algorithm is arranged from 
quasilinear QD chains considered in Section 2. We conclude our study by Section 4. 

2. The conditional electron transfer along a quasilinear chain of quantum 
dots 

In this Section we study a way of how to manipulate the two-electron system composed 
of charge qubit and auxiliary structure by selective driving of one of its components, 
namely, of the probe electron localized in the AS. As we shall see, the probe electron 
evolution in the AS depends on the qubit state, and, in its turn, may change internal 
qubit phase. Both aspects are important for an algorithm of entanglement production 
in many-qubit systems. 

2.1. The model and stationary eigenstates 

Consider the single-electron AS composed of electrostatically defined QDs {k = 
1 — A^) that are stacked in a quasilinear chain and separated by the finite potential 
barriers (figure 1). The structure parameters are chosen in such a way that the edge QD 
Ai {Ai\f) contains at least two bound electron states, namelly, the well-localized ground 
state \gi) {\gN)) and the excited state |ei) (Icat)) lying closely to the barrier's top, while 
each internal QD A^ [k = 2 — {N — 1)] contains at least one bound electron state |efc). 
The QDs Ai and A^ are supposed to be spatially isolated from each other so that the 
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electron tunneling between the states \gi) and Ig^) is negligibly small. If an electron 
occupies one of those states, it can stay there for extremely long time. We assume that 
the energies ^(ei) and ^(ejv) of excited states of the edge QDs and the energies e{ek) 
of states of the internal QDs are approximately equal to each other. Given the tunnel 
coupling Tk (Tfc > 0) between the states \ek) and \ek+i) is large as compared with the 
energy difference \e{ek) — e{ek+i)\ (resonant tunneling condition), those states become 
hybridized. If this condition holds for any neighboring AS QD states \ek) and \ek+i), 
then the single electron tunneling between QDs results in formation of a A^-fold excited 
state subband. Each state in that subband is a superposition of states \ek) {k — 1 — N) 
and is thus delocalized over the AS. Here we imply that the ground states of edge QDs 
do not hybridize with the excited states as long as £(ei(jv)) — e(5'i(Ar)) 3> ri(Ar_i). 

If charge qubit is positioned along the structure shown in figure 1, the 

energies of the AS QD levels are shifted relative to their unperturbed values s{gi), 
£{gN), and e{ek) by Ug{gi), Ug{gN), and Uq{ek), respectively, due to the electron- 
electron interaction. For definiteness, those shifts are completely associated with bare 
single-electron levels of the AS QDs. The shifts depend on the logical qubit state \q) 
{q — 0, 1) and result in the suppression of the single electron tunneling between the 
neighboring AS QDs and A^+i provided that the electrostatically induced energy 
mismatch for individual QD states is larger than the corresponding tunneling matrix 
element, i. e., \Uq{ek) — Uq{ek+i)\ > r^. In order to eliminate the electrostatic energy 
shifts Uq{ek) of the exited QD levels and to recover the resonant character of the 
electron tunneling along the AS, one can apply to each AS QD the voltages that 
generate the energy shifts SUq{ek) = —Uq{ek). Obviously, the sets {^^o(efc)}^i and 
{6Ui{ek)}k=i of compensating voltages for two logical charge qubit states |0) and |1) are 
different. For example, if one uses the voltage set corresponding to the qubit state "one" 
and \[Ui (ejfc) — Uq (cj^)] — [Ui (e^+i) — Uq {ek+i)]\ > Tk for some k, the unperturbed AS 
excited energy subband is reconstructed only if the qubit state is "one". Hereafter, 
we shall neglect the dependency of tunneling matrix elements on the qubit state 
assuming that corresponding variations in the interdot barrier heights have small effect 
on Tfe. In general, one can control the values of Tk with the voltages on the gates defining 
the barriers. 

We analyze our two-electron system within the tight-binding approximation and 
restrict ourselves to the case that the qubit electron is always localized in the logical 
subspace of the DQD. Thus the product states \gi,q) — \gi) \q), \gN,q) — Is'at) |?), 
ki, q) — |ei) ... , \eN, q) — |ejv) \q) {q — 0, 1) can be used as the basis states. In this 
case, the Hamiltonian describing the stationary two-electron states reads 

H — Has + Hq + Has-q, (1) 

TV 7V-1 

where i^AS = £{9i) \9i) (fl'i|+^(fl'7v) \9n) {9n\+J2 ^(efc) le^) {ek\- E \ek) {ek+i\ + h.c] 

k=l k=l 

is the AS Hamiltonian, Hg = |0) (0| |1) (1| is the qubit Hamiltonian, and the term 



HAS-q = 

9=0,1 



N 

Uq{gi) \gi, q) {gi, q\ + Ug(gN) \gN, q) {gN, q\ + Yl U'qi^k) \ek, q) (cfc, g| (2) 



fe=i 
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accounts for the effective interaction between the probe and qubit electrons. (In 
our model, the tunneling between the AS and the qubit as well as the non-diagonal 
electrostatic terms are completely ignored.) Here U'^^Ck) = Uq{ek) + SUq'{ek) and we 
assume that U'^{gi) ~ f/^(ei), Ug{giy) ~ Ug{eiy). The prime at U'^{ek) indicates on that 
the external voltages generate the shifts SUq'{ek) compensating the electrostatic shifts 
provided that q = q'. The electrostatic coupling energies are taken in the form 

^"""^ ^ l + {l-q)L/ro+{N-k)re/ro' 
where Uq = 1/ro, Tq is the distance between the center of the AS QD Aj^ and the 
nearest minimum of the DQD corresponding to the qubit state "one", L is the distance 
between the potential minima of the DQD, and Tc is the distance between the centers 
of neighboring AS QDs. For simplicity, the AS QDs arc supposed to be equally spaced 
from each other. In what follows, we shall consider the AS composed oi N — 20 QDs 
with uniform intcrdot tunnel coupling r at exact resonance when e{ei) = ... = s{ei\f). 
The dependencies of Uq{ek) on the AS QD positions are presented in figure 2. Hereafter, 
we shall work with the effective atomic units 1 au = Ry* = m* Ry / me£^ for the energy, 
1 au = = meSaB/m* for the length, and 1 au —h/Ry* for the time, where Ry is the 
Rydberg energy, is the Bohr radius, nie is the free electron mass, m* is the effective 
electron mass, e is the dielectric constant, and k is the Planck constant. Note that the 
parameters indicated in figure 2 and used throughout the paper correspond to realistic 
GaAs QD system (Ry* — 6 meV, a% — 10 nm). 

The voltage-controlled energy shifts are parameterized by the expression 5Uqi{ek) = 
—Uq'{ek) [6V/Uq], so that the external voltage compensates the electrostatic shift Uq'{ek) 
of the energy level in the AS QD Ak when 6V approaches Uq. Direct numerical 
diagonalization of the Hamiltonian, Eq. (1), results in two groups of two-electron 
eigenstates each corresponding to the qubit's localization in one of two logical states. 
The basis states l^^i, 0), \gi^, 0), \gi, 1), and \gN, 1) are the eigenstates of the Hamiltonian 
since we neglect the tunnel coupling between the AS ground states \gi) and \gN)- 
Their energies E{gi,0) = s (gi) + So + [/^(s'l), E(gN,0) = £ (gN) + £o + U(j(gN), 
E{gi,l) = £{gi) + £i + U[{gi), and E{gN,l) = £{9n) + £i + U[{gN) depend 
on the qubit state, the external voltage, and the structure geometry. The excited 
eigenstates of the probe electron may be represented as the normalized superpositions 

N 

1"^) Q) — ^ Cmkq kfe, where m — 1 - N and g = 0, 1. Since we shall be interested 

k=l 

in the PET between the AS ground states \gi) and \gN) localized in the AS QDs Ai 
and Am through the resonant temporal population of one of the excited states |m, g), 
it is important to know the superposition coefficients Crn,i,q and C.rn,N,q reflecting the 
weights of that transport state in both edge AS QDs. Figures 3 (a) and 3 (b) illustrate 
the dependencies of absolute values of those coefficients for the states from the central 
part of excited subband (8 < m < 13, A'" = 20) on the parameter 5V/Uq when the 
compensating voltages remove the electrostatic shifts corresponding to the qubit state 
"one" {q' — 1). The eigenenergies E{m,q) are shown in figure 3 (c). All structure 
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parameters are those given in Fig. 2. As it is observed from figures. 3 (a - c), the 
interaction of the AS with the qubit does not change considerably the resonant character 
of single-electron tunneling in the AS for both q = and q = 1 cases for a given set 
of parameters bringing about almost uniform shift of one excited subband relative to 
another. The weight coefficients converge to their unperturbed values for g = 1 at 
the point 6V = Uq where one has |Cm,i,i| = \Cm,N,i\ = |CAr_„+i,i,i| = \CN-m+i,N,i\- 
For g = a very close coefficient behavior is seen at 6V^OAUo, however, with some 
deviations from the symmetric picture demonstrated by the plots in the case of g = 1. 
We have revealed that the states lying in the middle of the excited subband conserve 
their hybridized structure in the range < 6V < 2Uo, whereas the edge subband states 
(not shown) are very sensitive to the external voltage. In particular, the tunneling 
collapse described above manifests itself in that the edge subband state with m = 1 (m 
= A^) splits off the subband and transforms at quite large positive (negative) values of 
6V {\6V\ > 3Uo) into the isolated state of the QD An with ICi^at,,]^ 1 {\CN,N,q\^ !)• 
The influence of the electrostatic interaction between the probe and qubit electrons on 
the AS spectral properties becomes stronger with the decrease in tq and/or with the 
increase in L. 

2.2. The PET and qubit phase shift in three-level approximation 

Next we demonstrate how to implement the PET between the ground states \gi) and 
Ig^) of the edge AS QDs Ai and Ajy given that the qubit is in one of its logical states, 
say, in the state Furthermore, we shall show that the PET from \gi) to Ig^) followed 
by the reversal PET from l^fjv) to \gi) can produce the phase shift between logical qubit 
states. 

Initially, the probe electron resides in the ground state \gi) of the QD Ai, whereas 
the state |g) = cq |0)-|-ci |1) of the qubit is arbitrary. Therefore, our two-electron system 
is characterized by the state vector 

1^,) = {Co exp [-^Eig^, 0)t] |0) + c, exp [-^Eig^, l)t] |1)} \g^) . (4) 

Note, that ?7o(ei)~f/i(ei) (see figure 2) and E{gi,0) — E{gi,l)^eo — ei, therefore, 
the phase difference between two state components in Eq. (4) is mostly defined by 
the difference between qubit logical state energies. When the compensating voltages 
are turned on, the excited two-electron eigenstates |m, 1) (m — 1 - A^) 
transform at 6V = Uq into the unperturbed eigenstates with high transport properties. 
According to Refs. [18] - [22], to attain the indirect resonant PET (now conditioned by 
the qubit state |1)) between the states \gi) and Ig^) localized in the QD Ai and A^, 
one has to irradiate the AS either by a pair of resonant laser pulses if s{gi) ^ s{gN) 
(asymmetrical case) or by single pulse if e{gi) ~ si^g^) (symmetrical case). Below 
we consider the asymmetrical case, where the ground-state energies of edge QDs 
substantially differ from each other. The pulse frequencies ujq and uoi match the resonant 
transition frequencies u {gi,r, 1) = E{r, 1)—E{gi, 1) and u {g^, r, 1) = E{r, l)—E{gN, 1) 
between the localized states \gi, 1) and \gM, 1) and some excited (transport) state |r, 1) 



Charge qubit entanglement in an array of quantum dots 



7 



[for concreteness, Uq = u {gi,r, 1) and ui = u {gN,f, !)]• Further, the pulse strengths 
Eq and El have to be chosen so that the absolute values of the coupling coefficients 
A(5fi,r, 1) = Eod{gi,r,l)/2 and X{gN,r,l) = Eid{gN,r, l)/2, where d{gi,r,l) and 
d{g]\[,r, 1) are the corresponding matrix elements of optical dipole transition, be equal 
to each other, i. e., \X {gi,r,l) \ = \X {gN,r,l) \ = A. It is also important that only 
the state |r, 1) from the AS excited subband is to be selectively populated during the 
pulse action. It means that the differences A(r, r ± 1, 1) = E{r, 1) — E{r ±1,1) between 
the energy of transport state |r, 1) and the energies of the nearest states |r ± 1, 1) must 
be large compared with the coupling coefficient A: |A(r, r ± 1,1) | ^ A. Besides, we 
should prevent the population of states from another excited subband corresponding to 
the qubit state |0), by careful choice of the structure and pulse parameters (see Sec. 
2.3). Finally, it is assumed that \£ (gi) — e {gN)\ ^ A so that each pulse drives only 
its own transition. (This requirement is relevant only in asymmetrical case.) As soon 
as all above conditions are satisfied, only three states, viz. \gi,l), \gN,l), and |r, 1), 
are optically active, and the effective three-level Hamiltonian describing the resonant 
transfer process under the rotating-wave approximation (A ^ uJo,u!i) takes the form 

Hrwa = A (5-1, r, 1) \gi, 1) (r, 1| + A {gN,r, 1) l^fA^, 1) (r, 1| + h.c. (5) 

In the rotating frame, the coherent evolution of the state vector |\I'a) = ai^i \gi,l) + 
0'N,i\gNA) + o,r,i\r,l)y spanned by those states, is governed by the non-stationary 
Schrodinger equation id\^A)/dt = HrwaI'^a)- Let the pulses be switched on at t 
= 0. For the initial condition |\['a(0)) = \gi,l), the solution of this equation is well 
known (see, e. g., Ref. pO]) and describes the three-level Rabi oscillations: 

ai^i = cos^ i^B.t) , ajy^i = — sin^ (^^i?,t) , fir,! = 1= sin (2i7/jt) , (6) 

V2 

where Qr = \/V2 is the Rabi frequency. Thus, the complete PET from the ground 
state of edge AS QD Ai to the ground state of edge AS QD takes place in times 
Tji = (7i"/2 + iTn)/Qii, where n = 0, 1, 2, ... In what follows, we shall consider the 
shortest time Tq = 7^/20,^ as the PET time. 

At the end of the optical transfer {t = Tq), the state vector of our system in the 
laboratory frame transforms into 

1^) = coexp HE(f7i,0)To] |0) \gi) - ciexp [-lE^g^, l)To] |1) \gN) • (7) 

After the time tq during which no pulses act on the AS, we drive the probe electron 
from the state I^^at) back to the state \gi) in the same manner and then turn off 
the compensating voltages. As a result, the component of the two-electron state 
corresponding to the qubit state |1) acquires an additional phase 9 = 5E{Tq + tq) 
[in comparison with the free evolution case of Eq. (4)], where 5E = E{gi\i,l) — 
E{gi, l)\sv=Uo — ^ (.9^) — £ (gi), so that for t > 2Tq + tq the state vector reads 

1^;) = {Co exp [-lEig,, 0)t] |0) + c, exp {-tO) exp [-tE{gi, l)t] |1)} \gi) .(8) 

For 6 = 7r(2m + 1), m = 0, ±1, ±2, ... the above operations amount to the phase shift 
of the qubit logical state |1) by tt in the qubit frame - i. e., |g) \gi) {Z \q)) \gi), where 
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Z = |0) (0| — |1) (1|. It is worth noting that given phase shift is defined by the difference 
between energies of AS ground states rather than qubit states. If e (g^) — e (gi) = 
(e. g., when QD Ai and An are identical), the phase difference between degenerate 
states \gi) and Ig^) is not accumulated and we need another approach to produce a 
phase shift. For example, by varying the compensating voltage on the QD A^ (after 
the PET l^fi) \g]y) now driven by single pulse have been completed), it is possible 

TO 

to generate required phase shift 9 = J 5e' (f) dt' by appropriate choice of the voltage- 



controlled energy shift 6e' {t') of the state \gj\f) and the variation time Tq. The time Tq 
is thus exploited as the additional independent parameter to control the qubit phase in 
both asymmetrical and symmetrical cases. 



2.3. Numerical simulations of the PET in 2(N +2)-level case 

The formulas (5) and (6) describe an idealized three-level evolution through which PET 
is realized with the probability pjv i(To) = |aiv i(To)p = 1. However, other excited states 
|m, 1) with m 7^ r (especially those nearest to the transport state) participate electron 
dynamics resulting in some decrease of pn,i{To). As it was shown in Ref. [23], the main 
reason for such behavior is non-correlated off-resonant excitations of those states during 
the pulse action. The PET probability pjv i(To) accounting for these processes is now 
given by the expression p4] 



mi (^o) ~ 1 - E /m, 

fm = [2A ((71, m, 1)/A (m, r, 1)]' sin^ [ttA (m, r, 1) /2v^A {g^, r, 1)] , 

where A(m, r, 1) = E{m, 1) — E{r, 1) and we require that \X{gi, rn, l)|~|A(5'Ar, m, 1)| for 
arbitrary m. Moreover, despite of inefficiency of the PET between ground states l^^i) 
and \g]\f) of the edge AS QDs for two-electron state component in Eq. (4) corresponding 
to the qubit state "zero", off-resonant excitations from the ground state |5'i,0) to the 
states |m, 0) (m = 1 - A^) may also happen. It may amount to noticeable reduction of 
the probability pi^iTo) of probe electron to stay in the state \gi, 0) at the pulse end. 

The PET optimization is thus considered as the search of the pulse and structure 
parameters for which bothp7v,i(^o) andpi^o(^o) are as large as possible while the transfer 
time To is rather short. In particular, as we have mentioned before, it requires the careful 
choice of transport state |r, 1) that has to possess large coupling coefficient X{gi,r,l) 
and be well separated from neighboring states. According to the results of the work |24j . 
in a quasi-one-dimensional structure the states belonging to the central part of excited 
subband are the best candidates for this purpose. If the QD number A^ is even, they are 
the states with indices r = N/2 or r = N/2 + 1, whereas for odd A^ this is the unpaired 
state with the index r = (A^ + l)/2. 

In what follows, we consider a general AS+qubit state in laboratory frame 

1^) = E c.U.e-'^^^^''^^' \g,,q) + a^.^e-^^^^-"?)* \gN,q) + E S„,,e-*^(™''')* |m,g)|(10) 

9=0,1 I m=l ) 
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as the superposition of 2{N+2) two-clcctron eigenstates at SV = Uq, q' = 1. Since 
we suppose that there are no transitions between logical qubit states during the PET 
in the AS, the subspaces (ifi-i, 0) , l^fAr, 0) , |m, 0)^^^^^! and \^\gi,l) ,\gN,l) 
corresponding to different logical qubit states, are decoupled from each other and can 
be treated separately. To study coherent probe electron evolution in the AS under the 
influence of a pair of resonant square laser pulses polarized along the AS axis x, the 
non-stationary Schrodinger equation 

ida.g/dt = A(g) [cos (ujot) + cos (ujit)] a^, q = 0, 1, (11) 

governing the probability amplitude vector = (ai^g, ajv^, ai^g, 0^?^^) , is solved 
numerically for the initial condition aq(0) = (1, 0, 0, 0)^. Here the matrix A(g) 
with the entries h.kk'{l) = 2A(A;, k', q) exp [—iA{k, k', q)t] characterizes the optical dipole 
coupling strength for each transition between the two-electron eigenstates from given 
subspace, and indices k and k' run over all of those states. The intersubband coupling 
coefficients between the ground states \gi, q) , \gN, q) and an arbitrary excited state |m, q) 
are calculated within tight-binding model as 

A {gi{N),m, g) = Sfieidd (^i(iv), m, g) /2, d (g^N),m, g) = Cm,i{N),qdo, (12) 

while the intrasubband coupling coefficients can be found from the expressions 

>'{9i{N),9i{N),q) =TSfieid{N -l)rj4, X{gi,gN,q) = 0, 

A (m, n, q) = -e^^m E Cl^^k^fuM,, \[N - l)/2 - A; + 1] r,/2, ^^^^ 
k 

where Sfieu = ^o*bEq/ Ry* is the field energy (actually, the dimensionless field strength), 
do = {gi \ — X |ei) is the matrix element of optical dipole transition between the states 
\gi) and |ei) of isolated QD Ai, and the origin is placed at the axis x in the center of the 
AS. Note, that in Eq. (11) we suppose the condition Eq = Ei, underlying the coupling 
coefficients symmetry \\{gi,m, 1)| = \\{gi^,m, 1)|, to be fulfilled. One can derive Eqs. 
(12) and (13) from the optical dipole matrix element's definition substituting in it the 
superpositional form for an excited hybridized state |m,g). We neglect in Eq. (12) 
the interdot optical dipole transitions setting (^fil — x |efc^i) ^{gN\ — x Ick^^) ~ and 
express in Eq. (13) the level shifts as {ek\ — x \ek) = — [{N — l)/2 — k + 1] and 
(cfel — X \ek'^k) ^ (the latter are much smaller than those where k = k' and have been 
discarded). Though the intrasubband terms have large coupling values in comparison 
with the intersubband ones, their contribution into the dynamics is minor since they 
oscillate at very low frequencies too far from resonance with the driving pulses. The 
dipole approximation used in Eq. (5) and Eq. (11) holds if the effective AS length 
Ic — {N — VjTc is much smaller than the radiation wavelength X^. In our model Ic ~ 
10~^ m, ~ 10~^ m and, therefore, one has lc<^\w 

Setting g = 1, we find from numerical solution of Eq. (11) the maximal PET 
probability max(pAr^i) = PAr,i(To) and corresponding transfer time To as functions of 
e field and then calculate the probability j9i,o(^o) = |ai,o(^o)P- In figure 4 we show 
the dependency of max(pjv,i) (filled circles represent numerical data and solid curve 
visualizes the approximation of Eq. (9)) as well as the dependencies of pi^iTo) vs the 
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field energy Efieid for three values of re,: pi_o for tq = 2, p[q for tq = 3, and p'(q for 
To = 4 (solid lines with open squares). Here do = 0.22 (1 au = ea^), r = 10 (the index 
of the transport state), = 15.083, Ui = 18.417 and the remaining parameters are the 
same as above. We see that numerical and analytical data for max{pi\f^i) correlate well 
with each other thus confirming our arguments used to justify the application of Eq. 
(9). Prom the other hand, it is rather difficult to construct a reliable approximation 
for the probability pi_o- If ^ Tq, the oscillations of this function can be satisfactory 
reproduced by the off- resonant solution for two-level system involving the state |g'i,0) 
and excited state |m, 0) closest to the state |r, 1). Otherwise, if L <t; tq, the populations 
of both components of a two-electron state demonstrate very close behavior because now 
it is hard to distinguish between the qubit charge states and the selectivity requirement 
is violated. We are interested, however, in the situation where L and ro are of the same 
order. 

Comparing both plots we define the optimal values of £/ieid as an abscises of the 
graph points where the probabilities simultaneously achieve their maxima. It is directly 
observed from figure 4 that it takes place in the interval around £/ieZd = 0.08 {Eo^i 
500 V/cm for GaAs) where the probabilities of interest are max(pjv,i) = 0.9989, pi^ 
= 0.996 at Efieid = 0.077 and max(pAr,i) = 0.9955, p'/o = 0.9974 at e field = 0.08. The 
plot illustrating the transfer time Tq vs the field energy Efieid (figure 5) helps us to find 
corresponding values of Tq: Tq = 846 and Tq = 833, respectively. For GaAs one has 1 au 
= 0.11 ps and the transfer time is about 85 ps (that is by an order of magnitude smaller 
than coherence time expected in GaAs nanostructures). Note, that for field energies 
Afield ^ 0.12 the PET probability substantially reduces. 

Let us summarize the results obtained in this Section. We have observed that 
conditional PET in a quasihnear QD chain between the edge QDs Ai and can 
be carried out with high accuracy in relatively short times. If the probe electron is 
then returned back to its initial state (viz., the ground state of the QD Ai), the phase 
operation is performed on the qubit positioned near QD Ajv- From the other hand, if 
our quasihnear chain is a part of complex AS (the case considered below) and the QD 
Aj\f belongs also to the adjacent AS part, the PET between those AS parts becomes 
conditioned by the state of the qubit attached to the QD An- Therefore, the probe 
electron can be transferred from the QD Ai via the AS to any other qubit provided that 
this first (control) qubit is in the state "one". Thus, it makes possible to implement the 
controUed-phase operation on given pair of (remote) qubits. 

In general, to perform a phase shift by an angle 9 on the target qubit provided that 
the control qubit is in the state "one", we have to implement conditional PET from 
initial QD of the structure to the QD positioned near control qubit (if and only if its 
state is "one"). Secondly, the conditional transfer from this near-to-control-qubit QD 
to the QD positioned near target qubit (again, if and only if its state is "one") should 
be carried out. Therefore, if the two-qubit state is |11), phase is accumulated in target 
qubit according to the scenario given in this Section. If the two-qubit state isn't |11), 
the probe electron doesn't go to target qubit and, therefore, phase isn't accumulated. It 
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is equivalent to the transition from two-qubit state \ab) to state [exp(i5)]'^^ |a6), where 
a,b= (0, 1). 

3. Nine-qubit Shor encoding and error syndrome measurement 

It is known that any reahstic quantum computational scheme should be protected from 
quantum errors caused by the environmental decoherence and the external control 
imperfections. For this purpose, one may use the quantum error correction codes 
converting individual qubit state into the specific entangled state of several qubits that, 
after decoding, restores initial single qubit state (see, e.g., the book of Nielsen and 
Chuang [1] for detail). In this Section we present the algorithm realizing the nine-qubit 
quantum error correction code developed by Shor [23]. 

3.1. Auxiliary structure for Shor encoding algorithm 

The encoding procedure requires us to arrange the qubits in such a way that the 
interactions between them can be switched on/off on demand in controllable manner. 
Making use of the results of Sec. 2, we propose the following model of auxiliary structure 
to mediate the interqubit coupling and to assist the entanglement generation needed for 
the implementation of Shor encoding (see figure 6). Such an AS incorporates three 
branches T, T', and T" each represented by a T-shaped QD structure. The QD Ai 
enters all branches as their common first QD and contains the probe electron in its 
ground state \gi) at the beginning of encoding procedure. Consider in detail the upper 
branch T that is attached to the qubit whose state we are going to encode. Its vertical 
part replicates the quasilinear AS examined in Sec. 2. The QD An positioned near the 
central (control) qubit qa is used not only for operations on that qubit but also for 
operations, conditioned by the state "one" of the qubit qa, on two side (target) qubits 
qs and qc placed near the QDs and Cn- Thus the QD An serves as the common first 
QD for two quasilinear substructures with the symmetry axis perpendicular to the axis 
X. The parameters of four QDs Ai, An, Bn, and Cn are chosen so that their ground 
state energies satisfy following relations: e{gi) e{gN) ^{9bn) e{gBj^) = ^{qcn)- 
Other two branches T' and T" are obtained from the branch T through in-plane clockwise 
rotations around the center of the QD Ai by the angles 27r/3 and 47r/3, respectively. We 
see that each branch mediates the interaction between the probe electron and a cluster 
of three spatially separated charge qubits. The tunnel couplings between the branches 
are controlled by the gates G, G', and G" . When a negative voltage is applied to one 
of those gates, the potential barrier separating the QD Ai from corresponding branch 
grows exponentially, and the electron tunneling across the barrier quickly falls. As the 
calculations show the decrease of tunneling matrix element r between two neighboring 
QDs by two orders of magnitude is enough to treat those QDs as isolated from each 
other. We shall consider each blockage gate to be either opened (no voltage) or closed 
(the voltage sufficient for tunneling's suppression is turned on). 
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The AS parameters are taken so that to neglect the differences between the energies 
of all possible charge configurations in the cluster. Those energy differences can be 
minimized by appropriate choice of the QD number and/or the interdot distances in 
each branch. The cluster qubits are thus considered relative each other regardless their 
spatially distributed structure and, consequently, without any conditionality between 
them that might be caused by various configurations of the electron positions in 
the DQD structures. The spatial addressability of cluster qubits (understood here 
as the possibility of the state rotation of a chosen cluster qubit (s) without affecting 
the neighboring ones) is achieved due to the frequency and polarization selectivity of 
resonant pulses. 



3.2. Nine-qubit encoding scheme 

Let all qubits except the first qubit qa be initialized. We shall work in the qubit reference 
frame so that the phase multipliers arising from the qubit state energy differences will be 
omitted. We use the notation \qA<iB<ic) = \'iA)\'iB)\<ic) {lA^qBAc = 0, 1) for the three- 
qubit basis state of the cluster coupled with the upper QD branch T. Similar notations 
for the three-qubit states of remaining two clusters attached to the QD branches T' and 
T" are supplied with prime and double prime, respectively. In order to transform the 
AS+qubits state = [cq |0) + Ci |00) |000)' |000)" into the encoded state 

l^/) = [co|05) + Ci|l5)]|^7i), (14) 



where IO5) (II5)) = 2-^/2 [|000) ± |111)] 



|000)"± |111)" 



is the Shor 



e qubit state |0) (|1)), we need 



|000)'± |111)' 

code state (the code word) ^23^ corresponding to the sing 
to complete the following set of operations on the qubits and the AS. 

According to standard decomposition of the encoding procedure into the set of 
single- and two-qubit operations [1], one initially should perform two controUed-NOT 
(CNOT) operations where qubit qA acts as control qubit whereas qubits and q'^ are 
used as target qubits. Given operations entangle different clusters. Next, the Hadamard 
rotations are performed on each of those qubits. Finally, two CNOT operations are to 
be organized inside each of three clusters. In this case, central cluster qubit qA {q'A, q'A) 
functions as the control qubit whereas corresponding side cluster qubits are target ones. 
Those operations entangle the qubits belonging to given cluster. What should we do to 
implement all of these steps using the AS shown in figure 6? 



3.2.1. Entanglement between the clusters We perform the CNOT operation on the 
qubits qA and in following manner. Firstly, we have to transfer the probe electron 
from the ground state \gi) of the central AS QD Ai to the ground state \gN) of the 
QD A]^ provided that the state of qubit qA is "one". For this purpose, we create 
the effective quantum channel connecting the QDs Ai and Aj^ via the probe electron 
tunneling along vertical part of the branch T. The application of compensating voltages 
to the QDs Afc (A; = 1 - N) amounts to the formation of the quasilinear QD structure 
with high transport properties, as it was described in Sec. 2.1. Note that apart from 
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the electrostatic interaction between the probe electron and the electron bound in the 
qubit qa, Eq- 3, the energy level shifts in the QDs are also affected by other qubits. 
It is essential, however, that required PET is to be conditioned by the state of the qubit 
qa only. To make transfer process insensitive to total charge state of remaining eight 
qubits, the maximal difference between interaction energies (here, the difference between 
interaction energies of the probe electron occupying the QD with those qubits all 
loaded either in the state "one" or "zero") should be much smaller than the coupling 
coefficient A. In this case, the optical excitation of AS transport state will not depend 
on charge configuration of the qubits qb, ■■■ , q'c- We have found that this condition is 
satisfied, e. g., for Efieid ~ 0.1 and Tc = 3 when the QD number in the chain A^, 
Bn{Cn) is greater than 18 - 20. The contributions originated from those interactions to 
the compensating voltages can be expressed as averaged interaction energies for given 
QD. 

Additionally, we close the gates G' and G" in order to interrupt the tunnel coupling 
between the QD Ai and the branches T' and T" . The blockage gates are used here to 
minimize the number of individual QD states that participate in the formation of excited 
subband. Such a reduction is accompanied with an increase in the spacings between 
nearest hybridized states and, therefore, it enhances the selectivity of the driving pulses. 
The numerical calculations illustrating the dependency of the probe electron energy 
spectrum on the external voltages indicate on the dissociation of total AS spectrum 
into the set of states delocalized over the quasilinear chain formed by QDs A^ [k = 1 
- N) and the states pertaining to other AS QDs. Residual tunnel coupling between QD 
A^ and its neighbors in substructures attached to the qubits qs and qc has little effect 
on transport properties of our quasilinear AS and may be completely ruled out by the 
use of blockage gates analogous to those surrounding the QD Ai. In what follows, we 
shall not discuss in detail the formation of a quasilinear transport channel supposing 
that this task can be solved in all relevant cases. 

As soon as the channel has been prepared, the PET can be attained according to 
the two-pulse resonant scheme since e{gi) ^ s{gN) (see Sec. 2.2). At the end of the 
pulse action one has 

1^) = |*o) + 01$ (To) |100) |000)' |000)" \gN) , (15) 

where the state vector component corresponding to the qubit state "zero" is denoted 
by l^o) = Co |000) |000)' |000)" and $ (Tq) is the phase multiplier. This form of the 
state vector establishes the conditional probe electron evolution needed to construct at 
following steps the controUed-phase operations on two pairs of qubits, (g^, ^^^d (g^i, 
g^), where plays role of control qubit. 

Secondly, the Hadamard rotation H {q'^) = 2'^^^ [|0) + |1)] (0| + 2'^/^ [|0) - |1)] (1| 
is implemented on the qubit g^: 

1^) = H (g;,) l^o) + c^<l> (t) 2-'^' |100) (|000)' + |100)') |000)" Ig^) • (16) 

(The single-qubit Hadamard gate can be accomplished via optical technique described 
in Ref. [20].) Further, we should perform the phase operation Z on the qubit g^ via the 
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PET between the ground states \giy) and \g'j^) of the AS QDs Aj^ and A'j^, conditioned 
by the qubit state Iqa) = |1) (namely, the controlled-phase operation), without affecting 
the state component H (g^) |\l/o)- To do this, the gates G and G' are opened and two 
laser pulses with equal strengths and frequencies, polarized along axes x and x', act upon 
the AS. Optically active part of the AS is now presented by two identical quasilinear 
QD chains with the edge QDs Aj^ and A'^. (Some specific features of excited state 
hybridization in the simplest planar artificial molecule composed of three disk QDs 
were analyzed in Ref. [19].) The pulse frequency matches the resonant frequency for 
optical transitions connecting the ground states \gj\f) and \g'j^) with the transport state. 
If the ground-state energies of QDs Ai and An (A'j^) are substantially different from 
each other, the pulses do not address the ground state \gi) of the central QD Ai. From 
other hand, the AS QDs A^ and A'^ are identical, therefore, required phase shift has to 
be achieved through the voltage sweeping on the QD A'j^ (see Sec. 2.2). After necessary 
voltage manipulations and probe electron's return into the ground state \gN), the state 
vector becomes 



\^) = H {q'^) iM/o) + ci$ (t) 2-1/^ |100) (|000)' - |100)') |000)" M • (17) 

As it is clearly seen, above operations are equivalent to CNOT operation on the pair of 
qubits qA and followed by Hadamard rotation of the target qubit g^. 
Next we repeat above procedure for the qubit q'^ obtaining 

|vi>) = E{q':,)}i{q'^) |vl>o)+Ci$(t) 2-^100) (|000)' - |100)') (|000)" - |100)") M 

and then transfer the probe electron from the state \gN) back to the state \gi). 
Performing the Hadamard rotation H (g^) on the qubit qA we arrive at the state vector 



where 



^) = [co ^o) + 01$ (Ti) ^i)J , (19) 

^o(i)^ = 2-i (|000) ± |100)) (|000)' ± |100)') (|000)" ± |100)"). Note, that at 
this stage we have realized on the qubits g^, g^, and g^ the three-qubit Shor encoding 
scheme that is able to correct trivial phase error. It is easy to write the expression for the 
time Ti required for above operations as Ti = GTo + STnad + '^To, where we approximately 
set all transfer times to be equal to Tq, Tnad is the time reserved for Hadamard rotation, 
and To is the voltage sweeping time. The phase 9(Ti) = [e{gN) — ^(5'i)](^i — Tq) 
accumulated during the process has to be equal to 27rn (n is integer) in order to cancel 
the multiplier ^(Ti) = exp [—i9{Ti)]. It may be attained through the appropriate choice 
of the time Ti . 



3.2.2. Entanglement inside the clusters Finishing part of encoding algorithm consists 
in the implementation of CNOT operations on the target cluster qubits g^ (g^, g^) and 
qc {q'c, Qc) with the control qubit qA {q'A, g^)- Since X = B.ZB., we perform initially 
the Hadamard rotation on one of target qubits, then carry out the phase operation 
Z, conditioned by the control qubit state "one", on that qubit and, as the final step, 
perform again the Hadamard rotation on the target qubit. This scheme is illustrated 
by evolution of three-qubit cluster state \qAQBQc)- The phase shift by vr of logical state 
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"one" of target qubit qb (qc) requires the conditional PET from the QD Ai to the QD 
A]\f and then from the QD A^r to the QD i^^r (Cat) and back. The main steps that 
bring about the CNOT implementation on the qubit qs are given below: 

[|000) ± 1100)] \g,) 1000) \g^) ± |100) |^;v) ^ 

|0) [|0) + |1)] |0) \g^) ± |0) [|0) + |1)] |0) M ^^^^-z;^--^- (20) 

|0) [|0) + |1)] |0) \g,) ± |0) [|0) - |1)] |0) 1^^) "^"^ 1000) \g,) ± |110) \g^) . 

[Here we omit in expressions of the state vector the phase multiplier and 
normalization coefficients. The qubit states of other two clusters remain unchanged 
and are not shown.] After the repeat of CNOT on the qubit qc and the return of probe 
electron into the state \gi), the components of the nine-qubit state, Eq. (19), transforms 
into 

^o(i)) ^ 2-t (|000) ± |111)) (|000)' ± |100)') (|000)" ± |100)") . (21) 

This sequence of operations takes the time T2 = GTQ+iTHad+'^To- Having performed the 
set of CNOT operations on remaining two qubit clusters (g^, q'^, q'c) and (g^, g^, q'^) we 
arrive at the desired encoded state, Eq. (14). The total time Ttot of encoding algorithm's 
realization is now expressed by 

Ttot = Ti + 3T2 = 24To + IhTHad + 8ro. (22) 

3.2.3. Time and probability estimations of encoding procedure Let us evaluate the time 
Ttot using the value of Tq calculated in Sec. 2.3. Setting Tuad^'To^ 10 ps (the value that is 
longer than it was taken in Ref. [20]) we find from Eq. (22) Ttot^ 2.3 ns at Sfieid = 0.077. 
We see that the transfer processes are the most prolonged part of the Shor encoding 
algorithm. Hence, the probability of successful encoding implementation can be roughly 
estimated as pshor^'o^i^[{PN,i)'^'^, {pi,oY'^] given that all other steps were performed 
without errors. It yields us the value pshor^ 0.91 that is, of course, insufficient for 
large-scale fault-tolerant quantum computations. However, this estimation establishes 
the possibility to realize given scenario in proof-of-principle experiments. The way of 
further optimization of the proposed algorithm includes the search of pulse and structure 
parameters that would provide an increase of the probability pshor with simultaneous 
decrease of the time Ttot- 

3.3. Error syndrome measurement 

It is known that with the help of nine-qubit Shor code one may, in principle, correct 
arbitrary single-qubit quantum errors occurring independently as long as an appropriate 
error syndrome measurement strategy exists. Here we demonstrate an algorithm 
correcting the X-type quantum error, viz., the environmentally induced qubit state 
inversion. Suppose that such an error has occurred in one of cluster qubits g^, 9s, 
or qc- In order to detect and correct it, one should perform the measurements of the 
observables Z {qA) Z (g^) = [|00) (00| + |11) (11|] ® Ic - [|01) (01| + |10) (10|] O Ic and 
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Z (qa) Z (qc) = [|00) (00| + 111) (11|] 0Ib- [|01) (01| + |10) (10|] ® Ib (see Ref. [I] for 
detail). If the results of both measurements are "+1", the three-qubit cluster state 
hasn't been damaged. If the measurement of Z (qA) Z {qb) gives "+1" ("-1") and the 
measurement of Z {qAj Z {qc) gives "-1" ("+1"), the state of the qubit qc (qB) has 
been inverted. Finally, if both measurements amount to the value "-1", it indicates on 
the error in the qubit q^- Since the corrupted qubit is identified, one has to perform 
X rotation on that qubit to recover the encoded state [Eq. (14)]. The quantum phase 
error Z that is unitarily equivalent to the inversion error X can be corrected in the same 
manner. Note, that syndrome measurement enables us to define only the damaged qubit 
but not the qubit state, i. e., it conserves the entanglement. 

We propose the syndrome measurement technique that is based upon the same 
principle as the encoding procedure. Again, the probe electron evolution in the AS 
conditioned by the states of tested qubits helps us to reveal the quantum error. Consider 
the AS composed of two identical quasilinear QD chains and used to measure the 
observable Z (g^) 2' (g^) (figure 7). The central QD Aq belongs to both chains and 
contains the probe electron in its ground state Iqq) at the beginning of the syndrome 
measurement. Two edge QDs Aq and Aq with ground states Ig'^) and \gQ) are placed near 
the qubits and g^, respectively, and their single electron states are affected by qubit 
states. The QD parameters are chosen so that the energies e (go), e (g'o), and e {g'^) of QD 
ground states If^o), Wo)^ and Ig^) are equal to each other, i. e., e {go) = e {g'^) = e ((^g). 
The application of the set of voltages, compensating the electrostatic shifts of the AS 
QD energy levels when the qubits are in the state "one", results in three different types 
of energy spectrum. If both qubits are in the state "zero", an electron tunneling via the 
AS excited subband remains inefficient. If the state of one of the qubits is "one", the 
transport properties of the QD chain, attached to that qubit, are recovered. Finally, 
if both qubits are in the state "one", the excited states become delocalized over the 
whole AS. Next we irradiate the AS by laser pulse of the frequency that matches the 
resonant transition frequency for situation where only one qubit is in its state "one". 
The strength and polarization of the pulse have to be taken so that the dynamical 
properties of transitions [5^0) ^ k) ^ Id'o) and [5^0) ^ k) ^ Ifl'o) connecting the transport 
state |r) and the ground states localized in the central and edge QDs will be equivalent. 
It means that the transition times and PET probabilities have to be the same for two 
possible charge configurations corresponding to that case. As soon as those requirements 
are met, the pulse will produce the three-level resonant dynamics described by Eqs. (6) 
giving rise to the PET between the state \go) and only one of states Ig'o) and Ig'^) (of 
course, if only one of qubits is in the state "one"). Note, that in the case where both 
qubits are in the state "one", despite of high degree of hybridization of the AS excited 
states, the probability of probe electron excitation from the state l^^o) is very low because 
of the resonant transition frequency mismatch. As a result, the probe electron leaves the 
AS QD Aq at the end of pulse only if the state of tested pair of qubits is IqAQs) = |01) 
or \qAQB) = 1 10). Thus, an information about the eigenvalues of syndrome operator 
Z (g^) Z (gs) can be extracted from routine detection of the probe electron in the QD 
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Aq by measurement of current Iqpc via the quantum point contact QPC. The eigcnvahic 
+ 1 (-1) thus corresponds to the presence (absence) of the probe electron in the QD Aq. 
If the measurement indicates on that the probe electron has left the QD Aq, the pulse 
is applied again to return it into the state \go), and the compensating voltages are 
switched off. It is very important that we cannot distinguish between two equivalent 
paths (otherwise the projective measurement of qubit state would be realized). Prom 
other hand, if one arranges driving pulse so that those transitions become different 
enough (e. g., polarizing the pulse along one of the QD branches), the setup shown 
in figure 7 may be used for standard projective measurement of a charge qubit state. 
Other five syndrome operators can be measured in similar way. 

4. Conclusion 

In our paper we have proposed the way of the entanglement engineering in quantum- 
dot-based nanostructures. As an important practical application of the developed 
framework, we have considered (to our knowledge, for first time) the realization of 
the nine-qubit Shor encoding scheme in two-dimensional array of the semiconductor 
charge qubits. Highly-entangled Shor code states protecting the quantum information 
against an arbitrary single-qubit error are generated by means of appropriate optical 
and voltage manipulations on the combined system involving both qubits and the 
auxiliary structure. As it was shown, the PET (conditioned by qubit state) along the 
auxiliary structure together with single-qubit Hadamard rotations underly the encoding 
procedure. Numerical calculations confirm the possibility of successful implementation 
of the transfer protocol. In general, the indirect qubit-qubit coupling looks as the best 
solution to organize distributed quantum state manipulations in the charge qubit array. 
The described principle of entanglement production that makes use of an auxiliary 
system mediating the interactions among qubits can be also adopted for quantum 
information processing schemes with fully electrical control. 

5. References 

[1] Nielsen M A and Chuang I L 2000 Quantum Computation and Quantum Information (Cambrige: 

Canibrigc University Press) 
[2] Valiev K A 2005 Usp. Fiz. Nauk 175 3 

Valiev K A 2005 Phys. Usp. 48 1 (Engl. Transl.) 
[3] Oliver W D, Yamagnchi F, and Yamanioto Y 2002 Phys. Rev. Lett. 88 037901 
[4] Fabian J and Hohenester U 2005 Phys. Rev. B 72 201304(R) 
[5] Tsukanov A V 2005 Phys. Rev. A 72 022344 
[6] Contreras-Pulido L D and Rojas F 2008 Phys. Rev. A 77 032301 
[7] Qiao B, Ruda H E, and Wang J 2002 Journ. App. Phys. 91 2524 

[8] DiVincenzo D P, Bacon D, Kempe J, Burkard G, and Whaley K B 2000 Nature 408 339 

[9] Tanamoto T and Fujita S 2005 Phys. Rev. B 72 085335 
[10] Reina J H and Johnson N F 2000 Phys. Rev. A 63 012303 
[11] de Pasquale F, Giorgi G, and Paganelli S 2004 Phys. Rev. Lett. 93 120502 
[12] de Visser R L and Blaauboer M 2006 Phys. Rev. Lett. 96 246801 



Charge qubit entanglement in an array of quantum dots 



[13] Wcinstcin Y S, Hellberg C S, and Levy J 2005 Phys. Rev. A 72 020304 
[14] Tanamoto T 2000 Phys. Rev. A 61 022305 

[15] Schirmer S G, Oi D K L, and Greentree A D 2005 Phys. Rev. A 71 012325 
[16] Fedichkin L, Yanchenko M, and Valiev K A 2000 Nanotechnology 11 387 
[17] Tsukanov A V and Valiev K A 2007 Mikroelektronika 36 83 

Tsukanov A V and Valiev K A 2007 Russian Microelectronics 36 67 (Engl. Transl.) 
[18] Openov L A 1999 Phys. Rev. B 60 8798 

[19] Oh J H, Ahn D, and Hwang S W 2000 Phys. Rev. A 62 052306 
[20] Openov L A and Tsukanov A V 2004 Pis'ma Zh. Eksp. Tear. Fiz. 80 572 
Openov L A and Tsukanov A V 2004 JETP Lett. 80 503 (Engl. Transl.) 
[21] Paspalakis E, Kis Z, Voutsinas E, and Terzis A F 2004 Phys. Rev. B 69 155316 
[22] Tsukanov A V 2006 Phys. Rev. B 73 085308 
[23] Shor P W 1995 Phys. Rev. A 52 R2493 
[24] Tsukanov A V 2008 J. Phys.: Condens. Matter 20 315204 



Charge qubit entanglement in an array of quantum dots 
Figure captions 



19 



Figure 1. (color online) Auxiliary structure formed by quasilinear chain of QDs 
(circles) and supplied with controlling gates (rectangles) contains a single probe 
electron in the ground state \gi) of the QD Ai. The qubit (single-electron double 
quantum dot) is positioned at the right from the structure and loaded in the logical 
state "one". Below, the AS potential profile along the structure axis x is shown 
schematically for the cases where the compensating voltages are switched either on 
(solid line) or off (dashed line). 



Figure 2. (color online) Electrostatic shifts of the AS QD levels for two logical qubit 
states as a functions of QD index k (QD position). Here Ui(ei) — Uo{e\) = 0.00082 
au, Ui{eM) — Uo{eN) = 0.3 au, and ?7i(ejv) — C^i(ejv-i) = 0.3 au. 



Figure 3. (color online) Characteristics of the AS eigenspectrum vs the control voltage 

parameter SV/Uq. a) The absolute values of the weight coefficients Cm.i,q in the QD 
Ai. b) The absolute values of the of weight coefficients Cm,N,q in the QD Ajv. c) 
The excited AS eigenenergies E{m,q) shifted by £(ei) + Sq. Plots for q = 1 (0) are 
presented by thick solid (thin dashed) curves. The point SV = Uq is marked by the 
vertical dotted line. 



Figure 4. (color online) Populations max(pAr_i) andpi^o, Pi o, Pio vs the field energy 
s field- 



Figure 5. (color online) Transfer time Tq vs the field energy e field- Thin curve 
with filled circles shows the numerical data and solid curve visualizes the formula 
To = 7r/2f2H. 



Figure 6. (color online) Schematics of the auxiliary structure used in the Shor 
encoding implementation (see text). 



Figure 7. (color online) Auxiliary structure proposed for error syndrome measurement 
on pair of qubits qA and qs (see text). 



Charge qubit entanglement in an array of quantum dots 




Figure 2 



Charge qubit entanglement in an array of quantum dots 



o 




Figure 3 (a) 



0.35 



— 0.25 



o 




0.15 



0.05 



Figure 3 (b) 



Charge qubit entanglement in an array of quantum dots 



22 




Figure 4 



Charge qubit entanglement in an array of quantum dots 



23 




»<^"n A"k 



1 B 



G" 



>1 



An 



C'n 



Figure 6 



Charge qubit entanglement in an array of quantum dots 



Jqpc O 

ox 

O 

#o 

^<cxj> o o 

Figure 7 



